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Abstract: An excess of events with a rapidity gap between jets, over what would 
be expected from non-diffractive processes, has been observed at HERA. A process 
based on a perturbative QCD calculation of colour singlet exchange has been added 
to HERWIG. With this addition, HERWIG is able to describe the number of events 
with a gap between jets over the number without a gap. This gap fraction is 
predicted to rise at large rapidity intervals between jets which would only be visible 
if the detector coverage were increased. 




Figure 1: Feynman diagrams of processes involved, (a) and (b) are examples of the direct 
and resolved contributions in photoproduction. (c) is the colour singlet exchange, (d) is the 
representation of how an individual event would appear in rj-cf) space. 



At HERA the incoming electron is accompanied by a cloud of photons. The spectrum 
of photon virtuality, P 2 , is dominated by almost-real photons and, provided P 2 is prevented 
from being too large, the incoming photon can be regarded as real and the types of event it 
participates in are classified as photoproduction events. We use the usual experimental cut of 
P 2 < AGeV 2 . If the event contains one or more jets it is assumed to have a hard interaction 
and is then perturbatively calculable in QCD; hence it is referred to as a hard photoproduction 
event. Here we will refer only to these kinds of events. These can be of two main types: 
direct and resolved photoproduction, shown in figures 1(a) and (b). In the direct case the 
whole photon interacts with a parton in the proton. In the resolved case the photon 'resolves' 
into partons, one of which then interacts with a parton from the proton. If the transverse 
momentum exchanged is high enough, outgoing partons give rise to 'jets' of particles in the 
detector. 

In most cases of parton-parton scattering a coloured object is exchanged (e.g. fig. 1(b)) 
with an associated colour flow between the outgoing partons and the remnant particles. This 
can be modelled by a colour string which is stretched across the central rapidity interval. This 
string then fragments into particles which occupy the region between the two jets. However it 
is also possible for the exchanged object to be a colour singlet (fig. 1(c)). In this case the colour 
strings connect each outgoing parton with the remnant jet closest to it in rapidity. This leads 
to a suppression of particle production in the rapidity region between the two jets. It was first 
observed at pp colliders [|1| and subsequently measured at HERA ||. 

The presence of high Et jets in these events guarantees a high value of t and assures us of 
the applicability of perturbative QCD, — £/Aq CD ^> 1 (s is the centre of mass energy and t the 
invariant momentum transfer squared of the parton system). In addition, the pseudorapidity 
interval between the jets, Arj, reflects the separation in rapidity of the outgoing partons, y. 
Therefore at large At] we have y ~ ln(— s/t) and —s/t ^> 1. Gluon exchange in the t-channel 
increases with 1/t 2 . However because a gluon is a coloured object, a rapidity gap does not 
normally result. Two gluon exchange, which may be in a colour singlet state, can give rise to 
a rapidity gap. Gluons exchanged between this t-channel gluon pair can further enhance the 
cross section such that it rises faster than 1/t 2 at small t. These gluons can be summed using 
the BFKL equation 0, to give 



as derived in ref. || {Cp = 4/3, Ca = N c = 3, ((3) ~ 1.202 is the Riemann zeta function and 
— c '^ as 41n2). The approximations made in deriving this formula mean that the correct 
scale to use in a s cannot be determined, and does not even need to be the same in each case. 
We use a s (—t) in the prefactor, a s = 0.25 in the denominator and u = 0.3 (which have been 
installed as HERWIG defaults). 

The package HzTool was used to generate the Monte Carlo and to compare to the 
original analysis || on ZEUS 1994 data. This facilitated the comparison of the data to the 
Monte Carlo generated. HERWIG 5.8d M was upgraded to include the QCD calculation |J 
described above. About 2.6pb _1 (the 1994 ZEUS luminosity) of events were then generated. 
A cone-based jet finding algorithm was run on this sample. Events were required to have at 
least two jets of E T et > 6 GeV with pseudorapidity satisfying 7]j et < 2.5. Denoting by r/i and r) 2 
the pseudorapidities of the two highest Et jets, the events were required to satisfy in addition 
| 'T 1 ^ 2 1 < 0.75 and At] = rji — i] 2 > 2. Events with no particle with transverse energy greater 
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than 300 MeV between the two jets were then classified as 'gap' events. The characteristic 
signature of these events is illustrated in fig. 1(d). 

The selected sample exhibits an exponential decay in transverse momentum (fig. 2(a)) and 
a bias towards high x° BS (fig. 2(b)). x® BS is the fraction of momentum of the photon that 
participates in the interaction as calculated from the observed jets. The differential cross section 
as a function of Arj is also shown (fig. 2(c)). The gap fraction (fig. 2(d)) is the number of events 
with a gap divided by the total number of dijet events. It levels off at around 60%. Although a 
fraction of 100% might be expected, many gap events are lost by such factors as the final state 
particles escaping the bounds defined by the jet cone and hence filling the gap. 
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Figure 2: The QCD colour singlet process by itself, (a) shows the spectrum of the jet used 
in the sample, (b) shows the x° BS distribution of the sample, (c) shows the cross section of 
the sample as producd by HERWIG, in bins of Ar] (the interval between jets), (d) shows the 
fraction of events with a gap over the total number of dijet events. 



HERWIG 5.8d was then used to produce 'standard' direct and resolved photoproduction 
events and these were added to the sample of 'colour singlet' events. The cross section for 
the dijet sample without the requirement for a jet was tuned to the ZEUS 1994 data. This 
was then held fixed as the 'colour singlet' sample normalization was adjusted to fit the gap 
fraction graph. An overall normalization factor of 30 for the 'colour singlet' sample was found 
necessary to describe the data. This factor is allowed due to theoretical uncertainties in the 
value of a s as mentioned above. PYTHIA 5.7 [[Jj was then used to produce 'standard' events, 
as a comparison. 

Both Monte Carlos are compared to the data (fig. 3(a)). As can be seen from the last two 
bins, PYTHIA, which does not include colour singlet exchange fails to describe the data. How- 
ever HERWIG, with the new QCD process, provides a good description. Also included is the 



standard HERWIG plot (HERWIG 5.8d) without the colour singlet sample. This emphasizes 
the dramatic effect, at high At] of the new process on the gap fraction. 
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Figure 3: Gap fractions as produced by PYTHIA and HERWIG (with the colour singlet inter- 
action denoted by 5.8d+). (a) is for a standard detector, and the ZEUS 1994 data has been 
superimposed, (b) is for an extended coverage detector. 

Using this model the Monte Carlo was used to simulate an extended detector. The extended 
detector coverage was simulated by changing the r]j et cut to r)j et < 3.5. With the extended 
detector clear evidence for the colour singlet exchange for Arj > 4.0 is seen (fig. 3(b)). Figure 3 
also shows a characteristic rise in the gap fraction. The gap fraction of the process by itself 
is 60%. The diffractive cross section falls less rapidly as Ar] increases than does the 'background' 
non-diffractive cross section. Thus at large Arj we expect the gap fraction to rise towards 60%. 
Hence given an increased detector acceptance, we can see a clear signal for this kind of event. 

The extended HERWIG was also used to simulate increased luminosity with the current 
detector coverage. Although increased statistics would allow different bins of Et of the jets 
and of x° BS no such striking signal of the colour singlet exchange would be seen. So we conclude 
in favour of the extended detector coverage. 

Thanks to Jeff Forshaw for useful discussions. 
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